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Abstract 



A generalized inverse problem for a two-dimensional difference opera- 
tor is introduced. A new construction of the algebro-geometric difference 
operators of two types first considered by I.M.Krichever and S.P.Novikov 
is proposed. 

1 Introduction 

The notion of a finite-gap with respect to one energy level Schrodinger 
operator was introduced by B. A. Dubrovin, I. M. Krichever and 
S. P. Novikov in the paper |IJ. The paper of S. P. Novikov and 
A. P. Veselov || deals with some class of two-dimensional Schrodinger 
operators called potential operators. In this paper S. P. Novikov and 
A. P. Veselov solved the inverse scattering problem. In the paper |J, 
I. M. Krichever has introduced a similar theory for difference operators. 
The recent papers |B|, ^[ deal with different natural generalizations of 
two-dimensional difference operators defined on regular graphs and lat- 
tices. In particular in the paper jjj (see also appendix I in the paper 
JBJ) in the context of the discrete Laplace transformations S. P. Novikov 
introduced an important class of the difference operators on equilateral 
triangular lattice. These papers stimulated new research in this area (see 
the review pi). 

In the present paper we propose a generalized inverse problem and a 
new construction of two-dimensional algebro-geometric operators both in 
Krichever's and Novikov's classes. 

Let L be a two-dimensional difference operator (of order 2K) 
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(1) 



i , j , 1 1 1 < A" , I j I < K 



with periodic coefficients 



a n + N,m ~ a n,m+M ~ a nm- 

Consider a space of Floquet functions 

V'n+iV.m = Wilpn.m, i>n,m+M = W2lpn,m- 

This space is finite-dimensional and the operator L induces in this space 
a linear operator L(tui, w^). The characteristic equation of this operator 

Q{wi,w 2 ,E) = det(E- Id - L(w 1 ,w 2 )) = 

defines a two-dimensional algebraic variety M 2 . A point of M 2 corre- 
sponds to a unique eigenvector tpnm of the operator L 

(Llp) nm = Elp nm 

such that 4>oq = 1. All other components il> nm are meromorphic functions 
on M 2 . Consider a curve F C M 2 corresponding to the "zero-energy level" 

T = {w 1 ,W2\Q{w 1 ,w 2 ,Q) = 0}. 

The functions ip nm are meromorphic on V. 

We can consider the two following problems. 

1) The direct spectral problem. Find explicitly the "spectral data" of the 
operator L (i.e. a set of geometric data like a curve F, divisors of poles of 
tpnm etc.) which determines the operator L uniquely. 

2) The inverse spectral problem. Find explicitly the operator L using the 
"spectral data". 

Both problems are complicated. It is nearly imposible to solve either 
of them in a general case. We can, however, consider a generalized inverse 
problem which consists of finding a set of geometric data with the following 
properties: 

1) The set of geometric data defines uniquely a family of functions ^„ m 
defined on an algebraic complex curve F, 

2) These functions satisfy the equation Lip = for some operator L of the 
form (|), 

3) The operator L is uniquely defined by the equation Lip = and the 
coefficients oJf m can be found explicitly. 

This problem is solved for some particular operators in the paper ||. 
I. M. Krichever calls such operators "integrable" but we will use the term 
"algebro-geometric" . 

Our goal is to find algebro-geometric operators. We found two exam- 
ples which can be of interest. 

The first example is provided by operators of the form 

(Ljp^nm — CLnm 1pn — l,m "F ^nm 1pn-\-\ y m + C-nmyJn ,771 — 1 

+ d 

(2) 

A value of (Li/j) nm depends only on values of ip at the points 
(n — 1, m), (n + 1, m), (n, m — 1), (n, m + 1), (n, m) 
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which form a cross in the plane (n, m). We will call such an operator 
"cross-shaped". These operators were considered by I. M. Krichever in |J. 
Algebro-geometric operators of the form (Q) found by I. M. Krichever 
correspond to a curve F' C M 2 whose image under the projection on the 
E- plane is the whole _B-plane. The corresponding problem is Lip — Eip, 
where both E and ip are functions defined on F'. In the present paper we 
deal with a different type of algebro-geometric operators of the form (Q) 
which corresponds to the "zero energy level" curve. The corresponding 
problem is Lip = 0. 

The other example is more complicated and maybe more interesting. 
Consider a triangular lattice in a plane. We will use as coordinates triples 
of integers k, I, m such that k + l + m — 0. On such a lattice we can consider 
an operator of the form 

(Llp)klm = aklm1pk,l + l,m-l + bklm 1pk,l- 1 ,m + l + Cklm'>Pk + l,l-l,m + (3) 
+4lm^fc-l,! + l,m + fklm1pk + l,l,m-l + Qklm V'fc- 1 , 1 ,m+l • 

A value of (Ltp)kim depends only on values of ip at the points 

(k,l + l,m - 1), (k, I - l,m + l), (k + l,l - l,m), (k- 1,1 + l,m), 

(k+l,l,m- 1), (k - 1, 1, m + 1) 

which form a hexagon in the plane k,l,m. We will call such an operator 
"hexagonal". In this case our lattice is not rectangular, nevertheless we 
can consider the generalized inverse problem and solve it. 

As we have already mentioned, this class of the operators has been 
introduced in the context of the discrete Laplace transformation by 
S. P. Novikov @, H||. 

It should be remarked that our formulas in the sections ^, ^ are not 
unique. We can choose other singularity structures for ^-function (for 
example using a rotation of the plane (n, m) by ^ in the case of operators 
of the form (13)) and obtain other algebro-geometric operators. 



2 Notation and conventions 

We use the notations and conventions of paper j5J. In particular, our 
conventions are the following. A basis of cycles oi . . . , a g , bi, ■ ■ ■ , b g is 
chosen in such a way that 

a,i o aj — bi o bj ■— 0, aiobj=Sij, i,j = l,...,g, 

where g is the genus of a non-singular curve V. A basis of holomorphic 
differentials wi, . . . , uj b is choosen in such a way that 

j> U} k = 2m8j k , j,k = l,...,g. 

a j 

We define the Jacobian J(T) as C 3 /{2mM + BN}, where M,N G Z 3 , B 
is a matrix of b-periods of u>i 

B jk = <* ujk, j,k = l,...,g. 

Jbj 
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We denote by Qpq the Abel differential of the third kind, i.e. a differential 
with unique poles at the points P and Q and residues +1 and —1 at these 
points respectively, we denote by Upq the vector of 6-periods of flpQ, we 
denote by K, the vector of the Riemann constants. 
We define the O-function as 

6(*)= E e*v(\{BN,N) + {N,z)Y 
New V ' 

where z = {z\, . . . , z g ) £ C 9 and ( , ) is a euclidean scalar product {x, y) = 

We use the following natural convention: if n is a negative integer then 
a zero (pole) of the nth order is a pole (zero) of the |n|th order. 



3 The cross-shaped operators: Krichever's 
class 

Consider an arbitrary two-dimensional difference operator L of the 
form (jij). Our goal is to find some solution of the generalized inverse 
problem stated in the introduction. 

Our construction is as follows. Let V be a nonsingular curve of genus 
g. Let P • , i = 1, 2, 3, be six points on V. Let D be a generic divisor of the 
form T> = Pi + . . . + P g such that the points Pk are different from the P*. 
Consider a function <f> a 0-y, ot, f3, 7 £ Z, defined on F such that: 

1) If a point P 6 V \ {P x , P^ 1 , P^} is a pole of (p a p-/, then P is one of the 
points P fc ; 

2) The function <^ Q( 3 7 has a zero of ath order in P x + and a pole of ath 
order in Pj~ , the same structure for /3 and P^ 1 , 7 and P^ . 

Lemma 1) Such a function a) 3 7 exists and is unique up to multiplication 
by a constant. 

2) The explicit formula for (j>a0~t i s: 

r Q/3r exp J^+Wi+lVs) Q {A (P)-A{V)-IC) ' 

where r af s 1 is an arbitrary constant, Po is a fixed point defining the Abel 



transform A (it should be remarked that the paths of integration in Jj 
and in the Abel transform are the same), fij = fl p + p - , Ui = U 



p 



p+p." 



Proof is done by a standard reasoning of the theory of the algebro- 
geometric integration. □ 

The key idea of the construction of our functions ip nm is a convinient 
relabelling in following way: ip n m = (frup-y, where 

a(n,m) = 2 ~ \~ m ,(3{n,m) = «-^™, 7 (n,m) = 
if n + m = (mod 2), 

a( n ,m) = 3 = \~ m , p(n, m) = ~ X +^ ~ m ,j(n,m) = 1+r ^ m , 
if n + m = 1 (mod 2). 

We will use a vectorial notation for the triples, i. e. the representation 
of a triple a, /3, 7 as a vector ai + f3j + 7k. For example, we will sometime 
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write a i+/3j+ 7 k instead of <j> a ,p,y This is useful because for exemple, if 
v = ai + /3j + 7k, then we can write 4> v +\ instead of <j> a +i,p,T 

We write ip nm = 4> v ( n ,m), where v(n, m) = a(n,rn)i + f3(n,m)j + 
7(n, m)k, i. e. 

v(n, m) = 2 ~" 2 ~ m i + + ^k, if n + m = (mod 2), 

v(n,ro) = 3 ~" 2 ~ m i+ -!+ 2 "-'" j + 1 +" 2 ~"' k, if n + m = l (mod 2). 

We will also use the following notation 

e(P,ai+/3j+7k) = 6(P,a,/3,7) = e(yl(P)+Q[7i+/3^2+7^3-^(2?)-/C). 

Let us formulate our theorem. 
Theorem 1. Let a family )/) ra „ be defined as stated above. Then Lip = 
if and only if the coefficients a nm , b nm , c nm , d nm , v nm of the operator L 
are defined up to a multiplication by a constant by the following formulas: 
1) if n + m = (mod 2), then 



a nm _ r v _j e(P+,v-j) . cxp( _ /0i) 



Pa 



d nm f v+i _j 9(P 2 + ,v + i-j) 

rv-j . e(P 2 +, v - j)0(Pf , v + i - j)0(P 3 - , v + i + k) 

dnm r v+k e(P 2 + ,V + i-j)e(Pf ,V + i + k)e(P 3 ",V + k) 

T P 2 + P l 

: exp( y Qi -J Qi (fi 2 + fi 3 )), 

Po -Po Po 

Cnm = r v _j e(P 2 + ,V-j)e(Pr,V + i-j) 

dnm r v+i+k 9(P 2 + , v + i - j)0(Pf, v + i + k) 

«nm _ r v _j e(P 2 + ,V-j)6(Pr,V + i-j 



X ( 

Po Po Po 

p i p r 

-Po Po 



dnm r v e(p 2 + ,v + i-j)e(pr,v + i + k) 

^ + p r p 2 

x exp(- / Qi -J (Sl 2 + Q 3 ) + J n 3 ) 



X 



e(p 3 -,v + i + k)e(p 2 -,v + k) 



n, r , e ( p 3 - ; v + k)e ( p.v) 

xexpc/no- ^r + V^ -cxpc/no' 
I e(Pa ' v) /o ■ 

where v = v(n,m), 

2) if n + rn = 1 (mod 2), then 

Anm = _ e(P 2 ~,V + j)9(P+, V - i+ j)9(P 3 +, V - i - k) 

Cnm r v _ k ' e(P 2 ",v-i+j)e(P 1 + ,v-i-k)e(P 3 + ,v-k) 

P 2 P l P 3 

x exp( j Qr+J (Q 2 + Sis) -J ^1), 
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Tv+j e(p 2 ,v+j) 



Cnm r v -l+j 9(P 2 )V — 1+j) 

dnm _ r v+j e(P 2 -,v+j)e(P 1 + ,v-i + j) 



exp( J Qi), 

exp ( y fii+y (n 3 +n 3 )), 

Fo Fo 



X 



e(p 3 + ,v-i~k)e(p 2 + ,v-k) . 
e(p 3 + ,v-k)e(p+,v) 



3 + 



c„ m r v _i_ k e(p 2 . v - i + j)6(P 1 + , v - i - k) 

Vnm = rv+j_ 9(P 2 ~,v+j)9(P 1 + ,v-i+j) 
Cnm r v e(P 2 ",v-i+j)e(P 1 + ,v-i-k) 

«T P l + P2 + 

x exp ( / fi a +J (n 2 + n 3 ) - y o 3 ) 

Po Po Po 

x exp(- / fix) - e ( p 2 + ,v -i-k) . ( _ / n j 

7 e(p+,v) PV 7 

Fo Po 
where v = v(n, m). 

Proof. Let Li/> = 0. Let us consider the case n + m = (mod 2). Thus, 
the formula Lip — becomes 

0V-J + "nm^v = 0. (4) 

Consider the point Pj~. Let A be a local parameter in a neighbour- 
hood of P]~. Hence a/ 3 7 = \~ a ■ h, where ft is a holomorphic function. 
The function exp fp Qi has a pole of first order at the point Pf . Thus 
exp Jp fl\ — K± A -1 + . . . , where if-j" is a constant. Therefore we have 



exp / f2 2 



V 



7 *r \ 7 



exp y n 3 



e(Pf, o,o,o) 



A- Q +. 



for P in the neighbourhood of Pj - . 

Now we can write down the term with / \-( Q ( n ' m )+ 1 ) m the formula 
in the neighbourhood of Pf" 



/ 



/3(n,m)-l 



7(n,m) 



, > a ( njm ) + l 
"nm 1 J ' v + i- 



V 



exp / f2 



exp y f2 3 



6(P 1 ,V + j-j) i-(ai( n ,n 

e(Pf, o,o,o) 



, ^ x - (a(n , m)+ D + Cnm( ^ i -)-(".-)+i rv+i+kX 




/3 (n,m) 



/ 



7(71, m) + l 



exp y f2 3 



9(Pf ,v + i + k) 

e(Pf, o,o,o) 



- (a(n,m) + l) 



= 0. 
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After simplifications we obtain a linear equation for the a nm and the Cnm 

a nm r v+ i_ J e(P 1 _ , v+i-j)+c„ m r v+i+k exp( J (Q,2+Q3))@(Pi , v+i+k) = 0. 

Pa 

Doing analogous computations at the points P 3 ~, P^ and P^ we obtain 
three other linear equations for the a nm , b nm , c nm , d nm and v nrn . These 
equations can be explicitely solved and the formulas for the coefficients of 
the operator L given in the statement of the theorem are obtained. The 
case n + m = 1 (mod 2) is analogous. 

Now let us suppose that the a nm , . . . , v„ m are given by the formulas 
of the statement of the theorem. Let us prove that Lip — 0. Consider the 
case n + m = (mod 2). Let us consider a function 

/ y i \ &nm I b nm Cnm , Vnm , 

(Llp) nm = Vn-l.m + ~j Wn + l,m + ~j Wn,m-1 + ~j Wn,m- 

dnm dnm dnm Qnm 

This function has the same pole or zero structure as the function —ip n ,m+i 
at the points P^ , i — 1, 2, 3. It follows from the formulas in the statement 
of the theorem that (Lip)„ m and — ip n ,m+i have the same pole or zero 
structure at the points P~ ,i = 1,2,3. If a point P G T \ {P*, P%, Pjt} 
is a pole of (Lip) nm or —tp n ,m+i, then P 6 V. Thus, by the Lemma, 
(Lip) nm and — t/j n ,m+i are proportional. Moreover, from the formulas for 
the coefficients of the operator L it follows that the terms with ^Z^"-" 1 ) -1 
in the series expansions of these two functions at the point P^~ are the 
same. Hence (Llp)nm ~ —tpn,m+i, but this is equivalent to Lip = 0. The 
case n + m = 1 (mod 2) is analogous. This completes the proof. □. 

Any set of non-zero constants g nm defines a "gauge" transformation 
of operators of the form (^) such that 

;_— 1 if — 1 i i -1 

a nm — 9n — l,m a ™ m 5 "nm — 9n+l,m^f lm ^ c nm — 9n,m — l c nm, 

,/ _ -1 , I _ -1 

This gauge transform acts on the eigenfunctions in the following manner: 

The following theorem is an easy corollary of the Theorem 1. 
Theorem 1' . For any set of "spectral data" consisting of: a non-singular 
curve L of genus g, six points P s 6 T, i = 1, 2, 3, and a generic divisor T> of 
g points different from the P*, there exists, up to a gauge transformation, 



a unique operator L of the form (^) 



4 The hexagonal operators: Novikov's 
class 

Consider a triangular lattice in a plane. We will use as coordinates triples 
of integers k, I, m such that k + I + m — 0. 

Consider an arbitrary two-dimensional difference operator L of the 
form Our goal is to find some solution of the generalized inverse 

problem stated in the introduction. 
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Our construction is as follows. Let V be a nonsingular curve of genus 
g. Let Qi, Ri, i = 1,2, 3, be six points on F. Let X? be a generic divisor of 
the form D — Pi + . . . + P g such that the points Pk are different from the 
Qi, Ri. Consider a function cj> a p-ypaT, a, /3, 7, p, a, r 6 Z, a + /3 + 7 = 0, 
p + a + t — 0, defined on T such that: 

1) If a point P € T \ {Qi, Q2, Q3, Ri, R2, R3} is a pole of a( 9 7po -T, then 
P is one of the points Pf, 

2) The function a /3 7 p CTT has a pole of ath order in Qi, a pole of /3th order 
in Q2, a pole of 7th order in Q 3 ; the same structure for p, a, r and Pi, 
R2, P 3 . 

Lemma 1) Such a function <^ a /3 7 p CTT exists and is unique up to multiplica- 
tion by a constant. 

2) The explicit formula for 4>ap-ypaT is: 

r a p~fpar ■ exp / (aQ,Q 3 Q 1 + I3Qq 3 q 2 + pQr 3 r 1 + <jQr 3 r 2 )x 
J Pa 

e{A{P) + aU Q3Ql + /3[/q 3 q 2 + pUr 3Ri + a?7fl3fl2 - ^(P) - g) 

6(A(P) - A(E>) - /C) 

where r a p JptTT is an arbitrary constant, Pq is a fixed point defining the 

Abel transform A (it should be remarked that the paths of integration in 

f„ and in the Abel transform are the same). 
Po 1 
Proof is done by a standard reasoning of the theory of the algebro- 

geometric integration. □ 

As in section |^ we will use a vectorial notation. We will represent the 

six integer numbers a, (3, 7, p, a, r as one vector 

v = aei + f3e 2 + 7e 3 + pe 4 + cre 5 + re 6 G Z 6 . 

Thus, we will write <f> v instead of (fiap-fpa-T- 

The key idea of the construction of our functions 4>kim is a convinient 
relabelling in following way: Vfcim = <t>v(k,l, m ), where 

v(fc, /, m) = — ei H — e2 H — e 3 H — — e4 H — e 5 H 3—66 

if fc-/ = (mod 3), 

v (fe, |, m ) = 4=i=lei + ^f±2e 2 + 2^±e 3 + *=|^e 4 + ^p±e 5 + 
+ m ' 3 fc ~ 1 e 6 if fc-Z = l (mod 3), 

v(k,l,m) = ^^-ei + ^f±±e 2 + 2^f^e 3 + £=^64 + ^f^e 5 + 
+ m-fc+i e6 jf fe _/ = 2 (mod 3). 

We will also use the following notation 

6(P, aei + /3e 2 + 7e 3 + pe 4 + ae 5 + re 6 ) = 6(P, a, /8, 7, p, a, r) = 

= 9(A(P) + aU Q3Ql + PUq 3 q 2 + pU R3Rl + oUr 3 r 2 - A{V) - K). 

Let us formulate our theorem. 
Theorem 2. Let a family ipklm be defined as stated above. Then Lip = 
if and only if the coeficients a k i m , bkim, Ckim, dkim, fklm, flfcim, of the 
operator L are defined up to a multiplication by a constant by the following 
formulas: 
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1) if fc Z = (mod 3), then 

Q>klm ^+64-65 



6(Q 2 , v - ei + e 2 )Q(Q 3 , v + e 4 - e 5 ) 
e(Q 2 , v + e 2 - e 3 )e(Q 3 , v - ei + e 2 ) X 



Qa Q2 

X ex P( J (^R 2 Ri - n QiQ 2 ) ~ J Q QsQi) + 
Po Pa 
| B(Q 2 , v + e 2 - e 3 + e 4 - e 6 )9(-Ri, v + e 4 - e 5 ) ^ 
6(Q 2 , v + e 2 - e 3 )6(i?i, v + e 2 - e 3 + e 4 - e 6 ) 



x exp( J Qr 3Ri - J {tt Q;i Q 2 + Qr 3 R 2 )) , 
Po Po J 

Qs 

dkim r v+e4 -e 5 6(Q 3 , v + e 4 - e 5 ) , [, 

Okim r v - ei+e2 B(y 3 ,v - ei +e 2 ) J 

Po 

fklm _ r v+ e 4 -e 5 8(-Rl, v + e 4 - gs) 



6tim r v+e2 _e 3+e4 _e 6 6 (Q 3 , v + e 2 — e 3 + e 4 - e 6 ) 
Ri 

xexp(/(-^ 3Q2 -^ 3fl2 )), £MHL =0 , ^ZH = o, 

where v = v(fc, Z, m), 
2)iffc-/ = l (mod 3), then 

&fcim _ 7~v-e 4 +e 6 9(-R 3 , V — e 4 + e 6 ) 



dfejm r v+ei _e 2 -e 5 +e 6 6 (i? 3 , v + ei — e 2 - e 5 + e 6 ) 
R 3 

X Cxp( j {£l Ql Q 2 + n Hl iJ 2 )), 



Cklm Tv — e 4 +eg 



dkim ^v+ei— e 2 



9(i?i, v + ei - e 2 - e 5 + e 6 )6(i? 3 , v - e 4 + e 6 ) 



9(i?i, v + ei — e 2 )6(i? 3 , v + ei — e 2 - e 5 + e 6 ) ' 

R 3 Rl 

x exp( J (£1 Qi q 2 + n_R lfl2 ) - J Qr 3 r 2 ) + 

Po Po 
| 0{Ri, v + ei - e 3 )9(Q 2 , v - e 4 + e 6 ) ^ 
9(i?i, v + ei - e 2 )9(Q 2 , v + ei - e 3 ) X 
Ri Qi 
xcxp( J £Iq 3 q 2 - J{ttQ 3Ql +Qr 3Ri )) 
Po Po 

Q2 

r v - e4+ee 9(Q 2 , v - e 4 



dkim r v+ei _e 3 6(Q 2 ,v + ei 

-Po 



~~~"T ' cxp( / (-fiQgQi - Jlflsfli)), 
- e 3 J J 
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0>kln 



gklr, 



dklm 

where v = v(fc, I, m), 
3)iffc-Z = 2 (mod 3), then 



bkl 



r v+e5 _ e6 9(Qi, v + e 5 - e 6 ) 



fklr. 



9(Qi, v - e 2 + e 3 ) 



exp( / (flQ 3 Q 2 + tt R3 R 2 )), 



Po 



e(i? 2 ,v + e 5 -e 6 ) 



J"v-e 1 +e 3 -e4+e 5 9(i? 2 , V — Bl + e 3 — e4 + e 5 



^2 

xexp( j \n QaQl +Qr 3Ri )), 
Po 



gklr. 



' v — e i +e 3 



0(03, v - e 2 + e 3 )9(Qi, v + e 5 - e 6 ) 
0(<?3, v — ei + e 3 )9(Qi, v - e 2 + e 3 ) 
Qi Q 3 
X exp(J(n Q3 Q 2 + Qr 3 r 2 ) - J tt QlQ2 ) + 
Po Po 
| 9(Q3, v - ei + e 3 - e 4 + e 5 )9(_R 2 , v + e 5 - eg) ;; 
9(03, v — ei + e 3 )0(i? 2 , v — ei + e3 — e4 + eg) 
R 2 Q 3 
xexp([(QQ 3 Q 1 +tl R3Rx ) + J fl Rl R 2 ) 



P 



O-klr. 



0. 



0. 



r. 

where v = v(fc, Z, m). 

Proof. The proof is analogous to the proof of the Theorem 1. □ 

Any set of non-zero constants hkim defines a "gauge" transformation 
of operators of the form (a) such that 



l t(m — h kl+lm _ 1 ahlm, b k i m — h kl _ lrn + 1 bl c lm, c klm — ^k+l,!-!, 



i (- k lr: 



dklm — h k _ ll + lm dklm, fklm = h k + llm _ 1 fklm, Qklm = h k _ llm + 1 gklm 

This gauge transform acts on the eigenfunctions in the following manner: 

V'fcim = hklmtpklm- 

The following theorem is an easy corollary of the Theorem 2. 
Theorem 2' . For any set of "spectral data" consisting of: a non-singular 
curve F of genus g, six points Qi,Ri 6 T, % = 1,2,3, and a generic divi- 
sor T> of g points different from the Qi,Ri, there exists, up to a gauge 
transformation, a unique operator L of the form (^). 



10 



Acknowledgements 



The authors are indebted to Professor Alexander P. Veselov for suggesting 
this problem and fruitful discussions. The authors also thank Professor 
Pavel Winternitz for discussions. The authors also thank K. Thomas for 
the help in the preparation of the manuscript. The authors wish to thank 
the referees for useful remarks. 

The main part of this research was performed during the participa- 
tion of one of the authors (A.O.) in the Seminaire de Mathematiques 
Superieures at the Universite de Montreal in the summer of 1999 and he 
is very grateful to the Universite de Montreal for hospitality. 

During this work the authors were supported by the grant INTAS 96- 
0770 (A.O.) and fellowships from the Institut de Sciences Mathematiques 
and the Universite de Montreal (A. P.), which are gratefully acknowledged. 

References 

[1] B. A. Dubrovin, I. M. Krichever, S. P. Novikov. The Scrodinger equa- 
tion in a periodic field and Riemann surfaces. Dokl. Akad. Nauk SSSR, 
229 (1976), no. 1., p. 15-18. (in Russian). English translation: Soviet 
Math. Dokl. 17 (1976), no. 4, p. 947-951. 

[2] A. P. Veselov, S. P. Novikov. Finite-zone, two-dimensional Schrodinger 
operators. Potential operators. Dokl. Akad. Nauk SSSR, 279 (1984), 
no. 4, p. 784-788. (in Russian). English translation: Soviet Math. 
Dokl., 30 (1984), no. 3, p. 705-708. 

[3] I. M. Krichever. Two-dimensional periodic difference operators and 
algebraic geometry. Dokl. Akad. Nauk SSSR, 285 (1985), no. 1, p. 31- 
36. (in Russian). English translation: Soviet Math. Dokl. 32 (1985), 
no. 3, p. 623-627. English translation: Russian Math. Surveys, 52 
(1997), no. 5, p. 1057-1116. 

[4] S. P. Novikov. Algebraic properties of two-dimensional difference oper- 
ators. Uspekhi. Mat. Nauk, 52 (1997), no. 1, p. 225-226. (in Russian). 
English translation: Russian Math. Surveys, 52 (1997), no. 1, p. 226- 
227. 

[5] S. P. Novikov, A. P. Veselov. Exactly solvable two-dimensional 
Schrodinger operators and Laplace transformations. In Solitons, ge- 
ometry, and topology: on the crossroad, Ed. V. M. Buchstaber and 
S. P. Novikov, AMS Transl. Ser. 2, 179 (1997), p. 109-132. 

[6] S. P. Novikov, I. A. Dynnikov. Discrete spectral symmetries of small- 
dimensional differential operators and difference operators on regu- 
lar lattices and two-dimensional manifolds. Uspekhi Mat. Nauk, 52 
(1997), no. 5, p. 175-234. (in Russian). English translation: Russian 
Math. Surveys, 52 (1997), p. 1057-1116. 

[7] B. A. Dubrovin. Theta-functions and nonlinear equations. Uspekhi. 
Mat. Nauk, 36 (1981), no. 2, p. 11-80. (in Russian). English transla- 
tion: Russian Math. Surveys, 36 (1981), no. 2, p. 11-92. 



11 



